Trace and antitrace maps for aperiodic sequences, their extensions and applications 
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We study aperiodic systems based on substitution rules by means of a transfer-matrix approach. In 
addition to the well-known trace map, we investigate the so-called "antitrace" map, which is the 
corresponding map for the difference of the off-diagonal elements of the 2x2 transfer matrix. The 
antitrace maps are obtained for various binary, ternary, and quaternary aperiodic sequences, such as 
the Fibonacci, Thue-Morse, period-doubling, Rudin-Shapiro sequences, and certain generalizations. 
For arbitrary substitution rules, we show that not only trace maps, but also antitrace maps exist. 
The dimension of the our antitrace map is r(r + l)/2, where r denotes the number of basic letters in 
the aperiodic sequence. Analogous maps for specific matrix elements of the transfer matrix can also 
be constructed, but the maps for the off-diagonal elements and for the difference of the diagonal 
elements coincide with the antitrace map. Thus, from the trace and antitrace map, we can determine 
any physical quantity related to the global transfer matrix of the system. As examples, we employ 
these dynamical maps to compute the transmission coefficients for optical multilayers, harmonic 
chains, and electronic systems. 

PACS numbers: 61.44.Br, 05.45.-a, 42.25.Dd, 71.23.Ft 



I. INTRODUCTION 

The trace-map technique, first introduced in 1983,0 
has proven to be a powerful tool to investigate the elec- 
tronic spectrum of various aperiodic systems, such as 
the Fibonacci sequence (FS)Jil the Thue-Mor.se sequence 
(TMS)J3 and the period-doubling sequence!] It has also 
been applied to investigate other nhysical systems, for 
instance, kicked two-level systems,Qtl and classical and 
quantum spin systems The technique was extended to 
study aperiodic systems in combination with the real- 
space renormalization-group technique. Recently, trace 
maps have been used to evaluate localization properties 
in a FS tight-binding modeLH 

This technique was transferred to the field of optics 
in order to see the scaling of the light transmission co- 
efficient through a Fibonacci dielectric multilayer .□ Now, 
we consider light that is vertically transmitted through 
a Fibonacci multilayer of two materials a and b which 
is sandwiched by two media of type a. The FS is con- 
structed by the substitution rule b — > a, a — > ab, The 
corresponding transfer matrices Ai are written as 





= PabPbPba, 




A 2 


= Pa, 




Ai+i 


= AlAl-i, 


(1.1 



where P a b{Pba) stands for the propagation matrix from 
layer a to b (b to a) and P a is the propagation matrix 



through the single layer a. They are given by 

Pab 



p-l 
r ba 



1 

n a /n b 



Pa = 



COs5 a 

sin <L 



- sin<5 a 
cos 5„ 



(1.2) 



where S a = kn a d a , n a is the refraction index of material 
a, d a denotes the thickness of the layers, and k is the 
wave number in vacuum. The quantity S a is the phase 
difference between the ends of a layer. For material b, 
the quantities 8b, rib, and db are defined analogously. 
The transmission coefficient is given byu 



\Ai 



(1.3) 



where \Ai\ 2 is the sum of squares of the four elements 
of Ai. Since the transfer matrix is unimodular, we can 
express the transmission coefficient as 

4 
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2 ' 



(1.4) 



where xi and yi denote the trace and antitrace of the 
transfer matrix Ai, respectively. Here, the so-called "an- 
titrace" of a 2 x 2 matrix 



.4 



An A12 
A21 A22 



(1.5) 



is defined as tja — A2 1 — Ayi, which follows the notion 



of Ref. 10. From Eq. (1.4), we see that the transmission 
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coefficient is completely determined by the trace and the 
antitrace, i.e., a complete description of the light trans- 
mission through general aperiodic multilayers requires 
both the trace and the antitrace map.E£l 

Now, we consider a different system, namely a har- 
monic chain composed of two kinds of masses, m a and 
TOfc, which are arranged according to the FS, and are cou- 
pled by two kinds of springs, K aa and K a {> = K\, a . Mak- 
ing use of tb p transfer-matrix formalism, the equation of 

iSS 



motion 



u n +i 



K 



Pa 



U n -1 



(1.6) 



where u n is the displacement of the nth atom from 
its equilibrium position, K n ^ n ±\ denotes the strength 
of the harmonic coupling between neighboring atoms, 



K 



+ 1 



m n uj 



and lu is the vibration 



frequency. 

From the corresponding global transfer matrix Ai = 
Y[ n =N t the transmission co effici ent ti and the Lya- 
punov exponent T/ are given b 



k 

Ti 



4 sin 2 k 



(zi cos k — yi) 2 + xf sin k 



yf - 2), 



(1.7) 
(1.8) 



where z\ = (Ai)n — (A;)22, Ni denotes the number of 
atoms in the chain, and cosfc — m a /(2K aa ). 

As our third example of physical systems based on 
aperiodic substitution sequences, we consider the trans- 
mission in electronic systems. This is closely related to 
the harmonic chain considered above. The Schrodinger 
equation for a one-dimensional tight-binding model with 
nearest-neicjhJpor hopping can be written in matrix form 
as followst3EJ 



= Mr, 



(1.9) 



where 4> n denotes the amplitude of the wave function in 
the Wannier representation, E the corresponding energy, 
e n the on-site energy at site n, and t rhn ±i the hopping 
matrix element between two neighboring sites. M n is 
the local transfer matrix associatedr-with site n. The 
transmission coefficient is given byt3t-a 



ti 



4-E 2 



(ziE/2- yi ) 2 + x 2 (l-E 2 /4)> 



(1.10) 



where the quantities xi, yi, and zi are again related to 
the global transfer matrix of the chain, i.e., the prod- 
uct Ai = n n =jv ^ n °f t ne local transfer matrices along 



the chain. Note the similarity to Eq. (1.7). The cor- 
responding Lyapunov exponent Ti is given by the same 
expression (] 



1. 



For the latter two systems, the Lyapunov exponent 
is completely determined by the trace and the antitrace; 
however, we need to know zi to calculate the transmission 
coefficient. Fortunately, it turns out that the maps for zi 
and yi are the same, as will be shown in Sec. IV. There- 
fore, the trace and antitrace map are sufficient to deter- 
mine the transmission coefficient and the Lyapunov ex- 
ponent. Thus it is desirable to construct antitrace maps 
for various aperiodic sequences, which is the motivation 
of the work presented here. 

The paper is organized as follows. In Sec. 0, we 
give the antitrace maps for various classes of aperiodic 
sequences, including the FS, the TMS, the periodic- 
doubling sequence, and certain generalizations. The ex- 
tension to arbitrary substitution rules and to maps for 
matrix elements are discussed in Sec. Ill and in Sec. IV, 
respectively. It is shown that the antitrace maps and 
the maps for matrix elements exist for arbitrary substi- 
tution rules and that the maps for non-diagonal elements 
and for the difference of the diagonal elements coincide 
with the antitrace maps. Applications to the computa- 
tion of transmission coefficients and Lyapunov exponents 
in different aperiodic systems are investigated in Sec. 0. 
Finally, in Sec. VI, we conclude. 



II. 



TRACE AND ANTITRACE MAPS FOR 
TWO-LETTER SEQUENCES 



We now proceed with the derivation of the antitrace 
maps of various classes of aperiodic sequences. We also 
include the corresponding trace maps for reasons which 
will become clear later. In this part, we make ample 
use of several relations for unimodular matrices. There- 
fore, we append a compilation of these relations in Ap- 
pendix ^. 



A. Generalized Fibonacci sequences 

There are many kinds of generalized FSs. Here, we 
study two- letter sequences F-SjCm, n) that can be gener- 
ated by the inflation schemeL-§EJ3 



<5o —b, Si — a, — St S) 



i 



(2.1) 



with arbitrary positive integers m and n, where FS(1, 1) 
corresponds to the well-known standard FS. Equiva- 
lently, they can also be generated by the substitution 
rule 

b^a, a^a m b n . (2.2) 

The total number of letters a and b in the word Si is 
denoted by Fi and satisfies the recursion relation 



F l+1 = mF x + nF Ml F = F 1 = l. 



(2.3) 
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In the limit of an infinite sequence, the ratio of word 
lengths for subsequent inflation steps is given by 



,. Fi+i 
nm 

I — >oo Fl 



+ V m 2 + 4n 



(2.4) 



Some values of a and commonly used terms for special 
cases of so-called "metallic means" ell are 



FS(1,1) : 




1 


f\/5 
2 


golden mean 


FS(2,1) : 


C s 


= 1 - 


I-V2 


silver mean, 


FS(3,1) : 




3- 


•h\/l3 
2 


bronze mean 


FS(1,2) : 


&c 


= 2 




copper mean 


FS(1,3) : 




1 


+ \/l3 
2 


nickel mean. 



It is known that the sequences FS(to, n) with n = 1 are 
quasiperiodic and those with n > 2 are always aperiodic. 

It is interesting to consider two further classes of gen- 
eralized FSsE 



~ x ab, 
- 2 aW 



The first class (2.5) a 



(2.5) 

2 ab. (2.6) 
ts of the so-called Fibonacci- 



class sequences FC(mji3o, the second (2.6) occurs in the 
renormalization-moup analysis of the energy spectrum 
of FC(m) chains.ta It is easy to check that the inflation 
schemes of these two generalized FSs are the same as 
those for FS(m, 1), but they differ in the initial words. A 
natural further generalization of these sequences is given 
by 



(b" 



(2.7) 



which we denote as FC(m, k). Her e, FC (m, 1) and 
FC(m, 2) correspond to the cases (2.5) and (2.6). The 
corresponding inflation scheme is 



Sq = b, Si 



&1+1 — <~>z 



(2. 



which is the same as that of FS(m, 1) apart from the 
different second initial word. 



1. The Fibonacci sequence 

Let us commence with the simplest example FS(1, 1). 
We consider the case that the two letters a and b corre- 
spond to two basic unimodular transfer matrices A and 
B, respectively. Denoting by A\ the total transfer ma- 
trix corresponding to a word S\ , the matrix equivalent of 
Eq. (O) for FS(1,1) is 



A t+1 = Ai-xAu 



(2.9) 



where A\ — A and Aq = B are the transfer matrices 
of the two building blocks a and b. Note the reversed 



order of matrix multiplicatio n as compared to the con- 
catenation of letters in Eq. (2.1), which occurs in the 
related tight-binding model that is usually considered, 
whereas the order of matrix multiplications is n ot re- 
versed in the optical problem, compare Eq. (jl . 1|) . The 
well-known trace map readsEJ 



Xl + l = Xl-lXl - Xi-2- 



(2.10) 



Note that in part of the literature a factor 1/2 is intro- 
duced in the definition of xi. Here, we omitted this fac- 
tor to keep symmetry between trace and antitrace. From 



Eq. (A.4), we obtain the antitrace map 
xiyi-i 



Vi+i 



yi-2- 



(2.11) 



The coefficients of the trace map are constants; however, 
those of the antitrace map include the traces. So, if we 
want to derive the antitrace map, the trace map must 
also be known. This is why we have to consider trace 
and antitrace maps at the same time. 



2. Generalized Fibonacci sequences FS(m, n) 



For FS(m,n), Eq. (2.1), the recursion relation for the 
transfer matrix is given by 

A l+1 = A\ 1 _ 1 A" 1 



(2.12) 



Here, we used Eq. ( ]A.l| ) and the definition of the func- 
tions U n {x) = C n -i(x/2) given in Appendix [a| in terms 
of the Chebyshev polynomials of the second kind C n (x). 
Furthermore, we introduced the notation 



U® = U n (x Al ). 



(2.13) 



From Eqs. fl2.12[ ), ( |A.3| ), and ( |A.4| ), the trace and the 
antitrace maps are obtained as 



Xl+l 



vi+i = 



yi+i = 



U'l + l = 



u {l - 1] u {l) v, -u (l ~Pu il) 

u n u m u l u n-l u m+1 



TT {l-l) T Al) 
" u n+l U m-1> 



r(0 



u n u m+l u l u n-l U m+2 

+ (x l+1 -UtPu ( J l l 1 )y l , 



(2.14) 
(2.15) 

(2.16) 

(2.17) 



where vi — XA t _ ± Ai and wi — ua,_,a,- No te th at the 
roles of vi and wi are subsidiary. Eqs. ( ^.14 ) and ( 2.15| ) 
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constitute the trace map, Eqs. (2.16) and (2.17) give the 
corresponding antitrace map. 

For special cases, these expressions simplify consider- 
ably. For FS(l,n), we obtain, using the properties (A. 3) 
of the functions U n (x), 



xi+i = U^vi-U^Zl ] x u 
v w = Ut^x m - uti } (*? 



-u 



(i-i) 

n+l I 



y l+1 = Ut^l-UtiVh 
wi+i = xi +1 y t + Uz^yi-x. 
Similarly, for FS(m, 1), we find 

xi+i = U$vi -U^^xi-i, 
vi+i = U^ +1 vi-U$xi-x, 



yi+i = u^wi-a 



lVi-i, 



Wl+l 



i rrO) rrO) 
xi +1 yi + U^wi - U^_ 2 yi-i 



(2.18) 

(2.19) 
(2.20) 
(2.21) 

(2.22) 
(2.23) 
(2.24) 
(2.25) 



Eqs. ( )2.22D -( |2.25| ) are quite different from Eqs. (pTTs|)- 
( 2.21 ) above. The corresponding aperiadic sequences 
show rather different physical properties .1121 We also point 
out that the trace and antit race m aps f or the sequences 
FC(m, k) are given by Eqs. (2.22)-(2.25) since they have 
the same inflation scheme. 

Eli minat in g th e subsidiary variables vi and wi in 
Eqs. ( 2.14 )— ( 2.17 ) for the general case FS(m,n), we ob- 
tain 



Xl+l 



yi+i = 



T Al) TT (l-l) 

Urn v n 



r('- 



u, 

T M) T Al-x) 

Urn 



1) {y-m+l x l u n+l u n-l J 



rjii) U {i-i)_ u (i) jtO-i) (2 26) 

u m+l u n-l u m-l u n+l ' l z - ZD ) 



LI 



(i-i) 



Ml) 7/0-1) „, . _ //(0 7/0- 1 ) „, 

u m+l u n yl-l u m u n-1 VI- 



(2.27) 

Here, we used Eq. ( |A.3| ) to simplify the result. The above 
two equations are alternative forms of the trace and the 
antitrace map. 

Again, for the special cases m = 1 or n = 1, these 
equations simplify. For FS(l,n), we find 



Xl+l — U n \U n _ 1 U n +1 J ' u n+l x l, 



(2.28) 



vi+i = ut 1} (ut 2) yi-2 + x m ^) -u«_p yi . 

(2.29) 

The result for FS(m, 1) reads 



U, 



0) 



Xl+l 



U, 



iy {u^+ixi - xi- 2 J - U^xi-i, 

(2.30) 



(') 



Um ( , rr0-l) N \ i rrO) 
Vl+l = \Vl-2 + U m-lVl ) + U^ +1 yi- 



(2.31) 



For FS(1, 1), Eqs. (p^)-(p^3l|) reduce to Eqs. ( pTl0| ) and 
( |2 . 1 1| ) , as expected. For some other special values of m 
and n, the trace and antitrace maps of FS(m, n) are given 
in Appendix Wi. 



B. Generalized Thue-Morse sequences 

Another type of aperiodiCj-aequence is the celebrated 
TMS and its generalizations. ErEil Here, we consider gen- 
eralized sequences TMS(m, n) with inflation schemeEU 

b^b m a", a^a n b m . (2.32) 

Equivalently, TMS(m,n) can be constructed as 

So = o, So = a, 



5Qm cm Q _ cm cm'. 

i+1 I / ) *->l+l — a l 



(2.33) 



For m = n = 1, this reduces to the standard TMS. 
The recursion relation for transfer matrices of TMS(ra, n) 
reads 



A l+1 = B?Ar, B l+1 = A\ n B l 



(2.34) 



where Ao is the matrix corresponding to the building 
block b, and Bq corresponds to a, respectively. 
Using the same method as above, we get 

arij.1 - U®tjW v , -7/W 7/0) _ 7/0) t/0) 

(2.35) 

- 770)7t0) 77O) 77O) 77O) 77O) 

- U 2n U 2m Vl - U 2n _ 1 U 2m+1 - ^2n+l ty 2m-H 

(2.36) 

where U; = xs t Ai- These two equations determine the 
trace map. 

It is somewhat more complicated to derive the anti- 
trace map because y^, 7^ ys r Wc d efine Vi = Va , and 
yi = y Bl . Then, from Eqs. (|2~3^) , (|Q|) , ( |A.12|) , and 
( A.13| ), we have 



yi+i 



yi+i 



Wl+l 



Wl + l 



(2.37) 

U$ {u^i-U^-U^U^yu 

(2.38) 

(u (l) U^v,-U m U {1) 

\ u 2n u m v l u 2n-l u m+l 

-UgUtU®^) U^yi + Ui%, (2.39) 

(n^U {l) vi-U {l) U (l) 

-U^Lult-i) U®v, + UW V1 . (2.40) 



4 



Here, Wi = VBiAi and Wi = i/a,b ,- T he anti trace map is 
completely determined by Eqs. ( 2.35 )-( 2.4p| ). 

For n — 1 and m = 1, Eqs. ( |2 .35[ )-( |2!40| ) reduce to 

(2.41) 
(2.42) 
(2.43) 
(2.44) 
(2.45) 
(2.46) 



xi+i = uj, 

= af (t7 f — 2) + 2, 

yi + i = w h 

m+i = m, 

= xi[(x l+1 - l)yi + yi], 

= xi[{xi+i - l)yi + yi]. 



m+i 

Wl + l 



This yields the well-known trace map of the TMS 

xi + i = x$_x(xi - 2) + 2, 
and the antitrace map 

vi+i = xi-x [{xi - l)yi-i + m-i] , 
m+i = xi-x [(xi - l)m-i + yi-i] ■ 

The above two equations give 

yi + i = xt-i[{xi +x/_ 2 - 2)y ; _i 

+X l - 3 Xi- 2 (2 - Xi- 2 )yi-z] 



(2.47) 

(2.48) 
(2.49) 

(2.50) 



which is an alt ernat iv e form of the antitrace map. 

From Eqs. ( 2.35 )-( 2.40 ), we can solve for the sub- 
sidiary quantities vi, wi, and wi, for instance, 



Vl+l 



u 2n u 2m 
u n u m 



xi + U. 



(') 



xi - U. 



(') 



-u. 



-u. 



-v. 



n-l \ u 2m x l u 2m-l J 
m-1 \ u 2n x l 



(I) 
n-l 

(0 



(2.51) 



The combination of Eqs. ( 2.35 ) and ( 2.51 ) gives an alter- 
native form of the trace map of TMS(m, n). 



C. Period-doubling sequence 

The period-doubling sequence can be generated by the 
substitution ruleu 



b — ► ba, a — ► b 2 , 
or the inflation scheme 

So = b, Si = ba, Si+i 



SrS, 



(2.52) 



(2.53) 



The inflation scheme is t he sam e as that of FS(1,2). 
Therefore, from Eqs. (2.26)— (2.27), the trace and the an- 
titrace maps are obtained as 



b l-2 



xi, 



xi +1 = jc/_i(x/afi_ 
yi+i = xi-iXi-^yi-2 + xm-iyi-i - y t . 



(2.54) 
(2.55) 



This yields also the trace and antitrace map for the cop- 
per mean sequence FS(1, 2). 



III. ARBITRARY SUBSTITUTION SEQUENCES 

As there-iestist trace maps for arbitrary substitution 
sequences, E3~Efl one natural question is whether antitrace 
maps also exist for arbitrary sequences. The answer is 
affirmative. We commence our argument in analogy with 
the discussion in Ref. |2^ and restrict ourselves to the case 
of unimodular matrices. 

Let At, A2, . . . , A r be 2x2 matrices and define the fol- 
lowing 2 r matrices 

B ei e 2 ...e r = A^A^ 2 . . . Ap (3-1) 

where ej <E {0, 1} for 1 < j < r. Then, from Eq. ( |A.S| ), 
any monomial Aj 1 Aj 2 . . . Aj s , with 1 < < r and 1 < 
i < s, can be written as_a linear combination of the 
matrices £? eie2 ... £r , namely,E3 



An Aj 2 



■ A, 



1 1 

y ] y ] ■ ■ ■ / , ' * i* ■■■■■■ ■ i-.- • 

£1=0 £2=0 e r =0 



c ^ 



(3.2) 



where each coefficient is a polynomial in the traces xa , 
1 < j < fi and the traces XA A k , 1 < j < k < r. 

This result not only yields the trace map, but also 
gives the antitrace map for any substitution sequence. 
We define 



Be 1 e 2 ...£ r ,l — A^A 



21 



(3.3) 



with ej € {0, 1}, 1 < j < r, and / > 0, where Aji is the 
unimodular 2x2 matrix associated to the l-th iterate of 
the j-th letter. Since each matrix in £? ei£2 ... ei .^ + i is, by 
definition, a monomial in the matrices A,;, they can be 
expanded in terms of the matrices i? eie2 ... er , j j according 
to ( |3.2| ). Then the trace of Z? eie2 ... Cr ,i + i is a polynomial 
in the 2 r — 1 traces of B 6l62 ,„ ert i; and the antitrace of 
B eie2 ... er .i+i is a polynomial in the 2 r — 1 antitraces of 
Be 1 e 2 ...e r ,b Therefore, we conclude that both the trace 
and antitrace maps exist for arbitrary substitution se- 
quences, and the dimension of the antitrace map is 2 r — 1 . 
Next, we present a concrete example to illustrate this 
conclusion. 

The Rudin-Shapiro sequence can be, defined by means 
of a substitution rule on four letters CJ The substitution 
rule is 

a — > ac, b — > dc, c — > ab, d — > db, (3-4) 

and the corresponding matrix recursion relations are 

A1+1 = C1A1, Bi+t = C1D1, 

Ci+t = BiA u D l+1 = B l D l . (3.5) 

We have the useful relationil 

Di = CiA^Bu (3.6) 

which effectively reduces the sequence to three basic let- 
ters. Now, we choose the seven matrices Ai, Bi, Ci, Di, 
A1C1, A1B1, and B\Ci as our basic set of matrices. 
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In what follows, we denote the traces and antitraces 



by 



a, 


— X A, 7 


b 




XB l , Q 


= x Cl 


j di = x Dl , 


ei 


= X A,C, 




fi 


= X Al B n 


9l = 


x B l C l , 


ai 


= VAo 


In 




VB X , Q 


= yc„ 


dl = VD t , 


h 


= UAiCi, 




f, 


= VAiBn 


m = 


y Bl c r (3.7) 



kDi 



By using Eqs. (AJ3) and ( A. 12 ), we obtain 
CiA u 

ciDi -aiBi + AiB h 
BiA h 

(aigi -cifi)Bi- ciAi 
-(gi-bici)AiBi 
+ {fi-aibi)QBi + C'iAi, 
fiCiAi - bid + CiB u 
c,[(l - af)Bi + eiDi + a t AiB t ] 
-CiB u 

hci[(fiai - afbi + b{)Ci 
+a l D l - C X B X - (/, - a,6,)C,A/] 
-h[{f l -aih){a l I-A l ) 
+kl + (af - l)Bi - aiMBi] 
-ciCi + I. (3.8) 

Note that the order of multiplication of two matrices on 
the right-hand sides of these equations may differ from 
the order of our basic ma trix products AiCi, AiBi, or 
BiCi. We can use Eq. ( A.S) to reverse the order to obtain 
a systems of equations that closes with our seven basic 
matrices. 



Ai+i 
Bi+i 
Ci+i 
Di+i 



Ai + iCi+i 
Ai+iBi + i 

Bi+iCi + i 



Now, from Eq. (3.5) and (A. 4), the trace and antitrace 
maps are obtained as 



bi+i 
Q+i 
di+i 
ei+i 
fi+i 

gi+i 

di+i 
bi+i 

Cl + l 

di+i 

ei+i 
fi+i 
9l+i 



c 4 di - aik + fi, 

fu 

bidi - ojcj + ei, 

eifi - b t ci + gi, 

ci(d t ei - afb t + aji + k) 

bici(aicifi - e t fi - afbia - 



-91, 

aibiei + kci 



+a t d t + gi) - 6, - c, + 2, 
cicii + aici - ei, 
cidi - aik + fi, 
bidi + aibi - fi, 

{aigi - aibid)bi + hdi - {gi - ha) fi 
+ (fi - aibi)(bici - gi) + a;Q - ej, 
fi{cidi + aici ~ e{) + c t bi - g t , 
ci(-afbi + edi + aji) - bih + g h 
6/(1 - cf)(fi - a t bi)ai + 6/(1 - of - cf)b t 



-cici + aikcidi + bici(fi - a ; 6/)e ; + a/6/// 
+6/c/gz. (3.9) 

Thus, we derived the trace and antitrace maps of the 
Rudin-Shapiro sequence. 

Now, we discuss the dimension of the antitrace map. 
Let A, B, and C be 2 x 2 matrices. ThenB 

ABC = [(xabc - xabxc - xaxbc + x A x B x c )I 
+ {xbc - x B x c )A - x A cB 
+ {xab - x A x B )C 

+x c AB + x B AC + x A BC}/2. (3.10) 

Taking the trace on both sides of Eq. fl3.10| ), we are led 
to a trivial identity. However, if we take the antitrace, 
we obtain 

VABC = [{x B C ~ X B X C )VA - XACVB 

+ (x AB - x A x B )yc 

+xcVab + xbVac + xaDbc] A (3.11) 

The antitrace of any monomial can be written as a lin- 
ear combination of a polynomial in the antitraces tja - , 
1 < j < f, and the antitraces yAjA k , 1 < j < k < r. Each 
coefficient is a polynomial in the traces xa 3 > 1 < j < f 
and the traces XAjA k , 1 < j < k < r. From this observa- 
tion we conclude that the dimension of our antitrace map 
is r(l + r)/2, i.e., the dimension is reduced from 2 r — 1. 
Here, for the dimension of the antitrace map, we do not 
take into account the dimension of the trace map, which 
enters the coefficients of the antitrace map. Thus, the 
full dimension of the trace and antitrace map is given by 
the sum of their respective dimensions. nn 
Let us consider two ternary sequences as examples E3'EII 
Our first example of a three-letter substitution rule and 
the corresponding recursion relation for the transfer ma- 
trices isB3 



a - 
A, 



6, 6 
i =Bi, 



c, 
B, 



■* ca, 

Ci, 



Ci. 



(3.12) 



Using Eqs. flA.lOj) and ( |A.12[ ), we obtain 

Bi, 
Q, 
Aid, 
CiBi, 

x Cl Aid - At, 
BiAid, 

BiAi + XBiAidCi 
-x BlAl I- 



Ai+i 
Bi+\ 
C/+i 
Bi+\Ai + i 
Ci+\Bi + i 
Ai+iCi+i 
Bi+iAi +1 Ci + i 



(3.13) 



Taking the trace of the above equation, we obtain the 
trace map. The dimension of the trace map is 2 3 — 1 = 7. 
Taking the antitrace of the sixth line o f the above equa- 
tion, we can expand it according to Eq. ( |3 .11 ). Therefore, 
the last equality in Eq. (3.13) is not necessary, and the 
dimension of the antitrace map is 3(3 + l)/2 = 6. 



G 



Our second example is the three-component FS gener- 
ated byEll 

a — ► b, b — ► c, c — > a&c, 

A m = B ; , B m = Q, C ;+ i = C1B1A1. (3.14) 
The corresponding maps for the matrices are 





-4/- 


-i 


= B h 






-i 


= C u 






-i 


= QBiAi, 






-i 


= CiB h 






-i 


= BiAi — XBiAj + Xc l BiA l Cl, 


Ci+iBi^ 




-i 


= Ai - xaJ + xctBtAiCiBi. (3.15) 



We see that, for this particular sequence, both the trace 
and antitrace maps are six-dimensional. 



IV. MAPS FOR MATRIX ELEMENTS 

As discussed in Sec. [| we need to know all elements 
of the global transfer matrix in order to compute certain 
physical quantities. Thus, the trace and antitrace maps 
may not be sufficient, and one would like to determine 
analogous maps for e ach of the matrix elements. 

Actually, from Eq. (3.2), we know that such matri x ele - 



ment maps exist for any substitution rule, and Eqs. ( 3. 
( |3.13 ), and ( 3.1 5| ) already contain examples of matrix ele 



ment maps. Now, we investigate the maps for the matrix 
elements of the FS (m ,n) an d TMS (m,ra ). 

Using Eqs. ([A|), flA.12| ), and ( |A.13j ), we obtain the 
matrix maps of FS(m,n) as (2.12) and 



+ (x l+1 -Uti 1) U«l 1 )A l 
-U r t^_ 2 A^ 

+ {utPu^_ 2 + v l+1 - x m+1 ) I, (4.1) 

where Mi = Ai^\Ai. The trac es xi +i and vi+\ appearing 
on the right-hand side of Eq. (4.1) are given in terms of 
Vi via Eqs. (2.14) and (2.15), respectively. 

Similarly, the matrix map of TMS(m,n) is obtained as 



A 



i+i 



UP U$N, 



V-i J 



-u. 



(') 



B, 



I m ' ^ m— 1 1 



-uWJuWAt-u^i 



r(') 



Ni+i = 



U (l) U^vi - U {1) U (l) 

u 2n u m u l u 2n-l u m+l 

~U {1) U {1) \ (u®Aj - U {1) , 

u 2n+l u m-l) \ u m A l u m-l J 



(4.2) 



(4.3) 



+ U 2n B l - U 2n+1^ 



(4.4) 



N, 



i+i 



in 



u n-l u 2m+l 



_Tj{l) TjU) \ 

u n+l u 2m-l J \ u n 



B, 



, TT (l) a TtW T 
+ U 2m A ! " U 2m+1 1 : 

where Ni = BiAi and Ni — AiB t . 



(4.5) 

We can eliminate the 



subs id iary matrices Mi, iVj, an d N i from Eqs. ( 2.12 ) and 
(4.1)— (4.5). For example, Eq. (4.1) becomes 



Mi = 



1 (tjV-V A, + IJV- 

Urn 

+VlI ~ XlXl-lI + XlAl-i. 



2) Ai-2 - Ut?l 



(4.6) 



Thus, we obtain anoth er fo rm of the matrix map of 
FS(m,7i) given by Eqs. fl2.12j ) and (O ). 

For m = n = 1, Eqs. §LA% and @ reduce to 



A 



i+i 



xixi-x)I + xiAt-x + A t -2- (4.7) 



This is the ma tri x m ap of the FS. For the TMS, we find 
from Eqs. (O)-(K^) for m = n = 1 



Ai+i = imKh-IJAhIBh-ih/]!/, 

(4.8) 

Bi+i = x l - 1 [(xi-l)Bi_ 1 +Ai_ 1 -xi- 1 I]+I. 

(4.9) 

The maps for the matrix elements are easily obtained 
from the matrix map, thus we do not give them exp licitly. 

Specifically, we consider the FS. From Eq. (p~7|), it is 
interesting to find that the maps for the non-diagonal 
elements, and for the difference of the d iagon al ele- 
ments, coincide with the antitrace map, Eq. (2.11). From 
Eqs. |[|) and (fD|, this fact also holds for the TMS. Ac- 
tually, as again follows from Eq. (3.2), we arrive at the 
important conclusion that the maps for the antitrace, the 
non-diagonal elements, and the difference of the diagonal 
elements are all the same for arbitrary substitution rules. 
This means that the knowledge of the trace and antitrace 
maps suffices to compute any physical quantities related 
to the global transfer matrix. 



V. APPLICATIONS 

We now turn our attention to applications of the dy- 
namical map method developed in this paper. In what 
follows, we are going to consider three examples. 



A. Optical Multilayers 

As our first example, we show how to use the antitrace 
map to calculate light transmission coefficients. 

The transmission of light through aperiodic multilayers 
arranged according to the FS,ErEa the "non-Fibonacci" 
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^equence,£^ the TMS,0 and the generalized TMSs0 
was studied in the literature. Possible applications of 
quasiperiodic multilayers as optical switches and memp* 
ries have been .suggested by Schwartz l3 Huang et ali3 
and Yang et aH found an interesting switch-like prop- 
erty in the light transmission through a FC(m) multi- 
layer. 

Using the antitrace map, we re-investigate the light 
transmission through FC(m) which is sandwiched by two 
media of type b. In analogy with the discussion of Ref. [n], 
we write the corresponding transfer matrices as 



A x = P b , 
A l+1 = A^A l 



(5.1) 



The recursion relation fo r th e transfer matrix (5.1) is a 
little different from Eq. (2.1) for FS(m,l). It can easily 
be seen that the trace map is the same, but that the 
antitrace map differs slightly. The antitrace map is given 
by 



yi+i 

Wl + l 



i rr(0 - rr(0 

xi + iyi + u^wi - u^_ 2 yi-i 



(5.2) 
(5.3) 



where w = y Al A l ^ 1 - 

We consider the case that the light vertically transmits 
the multilayer and choose the thicknesses of the layers 
d a and db appropriately in order to make fi a d a 
The n, w e have phase differences S a = S b = 5, compare 
Eq. (1.2). For S = (n + l/2)7r, the propagation matrices 
become 



Pa 



Ph 



-1 

1 



(5.4) 



Fro m the above equation and Eqs. (5.1), ( A.l), and 
flA.3p , we can obtain the initial conditions for the trace 
and antitrace maps as 



yi 

t/2 
W 2 



o, 

-U m -i(0)gi, 

U m - 2 (0)gi, 

2, 

-U m - 2 (0)gi, 
-J7 m _i(0)ei, 



where 



R m + R~ 



R = n a /n b . 



(5.5) 



(5.6) 



The initial conditions depend on the parameters Q\ and 
m, w hile the recursion relations only depend on m. From 
Eq. (A.3), we know that 



mo) = ^[i-c-in^- 1 





(-l) fc 



for m = 2k 
for m = 2k + 1 



(5.7) 



Since i7 m +4(0) = U m (0), i.e., the function U m (0) is peri- 
odic in m with period four, the initial conditions, given 
in Table Q. also show this periodicity. The initial condi- 
tions for FC(2g), q = 1, 2, 3, . . ., or for FC(2<? + 1), only 
differ by the sign of the parameter R. Thus, it is nat- 
ural to divide the FC(m) into two classes, FC(2q) and 
FC(2? + 1). 

From the initial conditions and recursion equations, 
we can directly obtain the trace, the antitrace and the 
transmission coefficients of FC(2g), which are given in 
Table |l| It can be seen that the trace and the antitrace 
vanish alternately. The trace shows periodicity with pe- 
riod four for odd values of q, and period two for even q, 
but the antitrace shows no periodicity. Thus, the trans- 
mission coefficient also is not periodic in /. For even the 
transmission coefficient does not depend on m. However, 
for odd I, the transmission coefficient depends on m and 
I, see Table |]. 

Table |l| shows the results for FC(2g+l). In this case, 
the trace, the antitrace and the transmission coefficient 
are periodic in I with period six. The transmission coeffi- 
cients are the same for I = 2, I = 3, and I — 6 and do not 
depend on m. We find that the multilayer is transparent 
for l = 6i + l, i= 0,1,2.... 

Here, we not only recover the recent results of Ref. [lj], 
but also give a natural classification of FC(m) and derive 
the periodicities of the trace and antitrace maps. 



B. Harmonic Chains 

As our second example, we show how to apply the 
map for the matrix elements to calculate some physical 
quantities for a harmonically coupled Fibonacci chain. 
The transmission coefficient an d the Lyap unov exponent 
were already given in Eqs. (1.7) and fll.Sj ). We know the 
trace ( 2.10Q and antitrace maps ( p. 11 ) for this system. In 
order to determine the transmission coefficient, we addi- 
tionally need to know the map for the difference zi of the 
diagonal elements in Eq. fll.7| ). As discussed in Sec. IV, 
the map for z\ is the same as the map for the antitrace 
Vi- 

Now, this leaves us with the problem to determine 
the initial conditipes. By a so-called transfer matrix 
"renormalization" the transfer matrix product can be 
re-written in terms of "renormalized" transfer matrices 
such that these are arranged according to the FS. Follow- 
ing the discussion in Ref. O, we choose a special value of 
parameters 



n = 



2(3+1 



a(l-/3) 



(5.8) 



where a = mb/m a and /3 = K a ^4,K ab . The first two 
renormalized transfer matrices aretil 



.4, 



1 

7/1 1 
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An 



-1 



>)2 



-1 



(5.9) 



where 771 = 2(a — 2) and 772 = 2(1 — a). Note that 
these two matrices commute with each other for arbitrary 
values of r\\ and 772. From this equation, we obtain 



A 3 = A X A 2 



-1 





-1 



Ji2-rji 

Thus, the initial conditions are given by 

xi = 2, x 2 = -2, x 3 = -2, 

2/1=571) V2=m, 2/3 = 772 — 771 , 

zi = z 2 = z 3 



(5.10) 



0. 



(5.11) 



From the antitrace map (2.11), we find that z\ = 
for all /. Using the trace map (2.1C), we easily obtain 
x 3l+ i = 2, x 3i+2 = -2, and x 3i = -2. That is, the 
trace map is pe riod ic in I with period three. Then, from 
Eqs. (1.7) and ([O]) the transmission coefficient and the 
Lyapunov exponent have the simple forms 



tr 



1 



vt 



4sin^fc 



T ; = TV -1 ln(yf + 2) 



(5.12) 
(5.13) 



From the initial conditions for y\ and the antitrace map 
( [2.11 ) , we easily find that the modulus of yi is 



(5.14) 



where Fi denotes the Fibonacci number defined by the 
recursion Fi = + F;_ 2 with Fq = F\ = 1. 

Finally, the transmission coefficient and the Lyapunov 
exponent are obtained as 



tr 1 = 1 



Ti = N 



4 sin 2 k 
- l \n[{F m - F^ im ) 2 



2]. 



(5.15) 
(5.16) 



Thus, using our matrix element maps, we have re-derived 
the result of Ref. 0. 



C. Electronic systems 

We now apply the trace and antitrace method to the 
transmission problem in electronic systems for the ex- 
amples of the FS and the TMS. In what follows, we 
choose the parameters as e a = — tb = e, tab = 1, and 

taa tbb — t. 



1. Fibonacci sequence 

For the FS, ther e ar e actually four different local trans- 
fer matrices M n ( |l.9| ), because the hopping matrix ele- 
ments depend on three subsequent letters in the FS. Nep-, 
crthcless, the transfer matrix product can be re-writtentil 



in terms of two matrices 



M a 



1 J \ 1 

E-e -t\ ( {E-e)/t -1/7 

1 oil 1 



E + e -1 
1 



(5.17) 



such that the resulting transfer matrix product is again 
arranged according to the Fibonacci sequence. 

For the trace and antitrace maps, we only need to know 
the first three matrices A\ = Ma , A 2 = MbM a , and 
A 3 = M a M b M a . From Eq. fl5.17 ), these matrices and 
thus the initial conditions are easily obtained. In order 
to obtain an analytical result, we restrict ourselves to the 
case E = e = 0. For this particular choice of parameters, 
Eq. fll.lOD simplifies to 



V! 



(5.18) 



which is formally the same as Eq. (1.4). The initial con- 
ditions become 



x\ =0, x 2 = 0, x 3 = 2, 

Vi = -t- 1/7, V2=t + 1/7, y 3 = 0. 



(5.19) 



From the trace and antitrace map equations (2.10)— ( p.ll ) 
for the FS, we can easily find that both the trace xi and 
the antitrace yi are periodic in I with period six. In one 
period the traces are 0, 0, 2, 0, 0, —2, and the antitraces 
are-t-l/f, 7+1/7, 0, t+l/t, t + l/t, 0. From Eq. (|5.18|), 



we deduce that the transmission coefficient 7; is periodic 
in I with period three. For one period, the transmission 
coefficients are given by 4/ (t+ 1/7) 2 , 4/(7+l/rj) 2 , and 1. If 
the hopping parameter 7=1, the transmission coefficient 
7; = 1 for all values of I, which is the trivial (periodic) 
case. Next we consider the electronic transmission for 
the TMS. 



2. Thue-Morse sequence 

We consider the on-site model for the TMS, i.e., the 
hopping parameter 7 = 1. So there are only two kinds of 
transfer matrices 



Bo 



E + e 
1 



-1 




A 



E-e 
1 



-1 




.(5.20) 



From these, we can calculate the matrices A±, B±, A 2 , 
and B 2 , and thus the initial conditions for the trace 
and antitrace map. Again, in order to obtain an an- 
alytical result, we limit ourselves to the case where the 
parameter e and the energy E fulfill a particular relation, 
E = y2 + e 2 . In this case, the initial conditions become 



•i'o 



= a/2 + e 2 - e, xi =0, x 2 = -2 - 4e 2 
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yo=yo = 2, V\=V\ = 2^2 + e 2 , 2/2 = -j/a = 4e, 

(5.21) 



-'o 



+ - e, 



z 1 = zi = 2, z 2 



"22 



wher e gj = (A ; ) n - (A;) 22 and z/ = (Bi) n -(Bi) 2T From 
Eq. ( p. 47] ), we deduce that the traces X[ = 2 for all / > 3. 
From the antitrace map equations ( 2.48 )-( 2.49 ) and the 
above initial conditions, we easily find that yi = zi = for 
/ > 3. Thus, we obtain the result that the transmission 
coefficient t; = 1 for I > 3. For I = 1 and 1 = 2, the 
transmission coefficients are given by t\ = (2— e 2 )/(2+e 2 ) 
and t 2 = (2 - e 2 )/(2 + 7e 2 + 4e 4 ), respectively. 

The examples considered here show that trace and an- 
titrace maps provide a convenient tool for the computa- 
tion of physical quantities related to the global transfer 
matrices of aperiodic substitution systems. In the ap- 
plications presented above, we mainly concentrated on 
obtaining analytical results, and therefore had to restrict 
the discussion to specific values of the parameters. The 
trace and antitrace map equations, of course, are not 
restricted to these cases, but there will be no simple 
closed-form solutions to the recursion relations in gen- 
eral. The particular parameter values considered above 
correspond to periodic orbits of the associated dynamical 
systems. These cases, and probably all examples where 
simple solutions exist, share the property that, at a cer- 
tain stage, different transfer matrices commute with each 
other, and thus are simultaneously diagonalizable. This 
also explains why these systems turn out to be transpar- 
ent, because it does not matter in which order one multi- 
plies matrices that commute with each other. In spite of 
these comments, the method presented here is expedient 
and useful for the investigation of physical systems built 
on aperiodic substitution sequences, because the trace 
and antitrace map equations can very efficiently be used 
in numerical investigations of large, but finite, systems. 



VI. CONCLUSIONS 

In conclusion, we have extended the well-studied trace- 
map method for the investigation of aperiodic substitu- 
tion systems by considering corresponding maps for the 
antitrace and the matrix elements of the transfer matri- 
ces. Our main results are the following. 

Firstly, we obtained the trace and antitrace maps 
for various aperiodic sequences, such as generalized FSs 
and TMSs, the periodic-doubling sequence, examples of 
ternary sequences, and the four-letter Rudin-Shapiro se- 
quence. The dimension of the dynamical systems defined 
by the trace map and our antitrace maps is r(r+l)/2 plus 
the dimension of the trace map itself, where r denotes the 
number of basic letters in the aperiodic sequence. Sec- 
ondly, we showed that trace and antitrace maps can be 
constructed for arbitrary substitution rules. Thirdly, we 
introduced analogous maps for specific matrix elements 
of the transfer matrix, but it turns out that the maps for 



the off-diagonal elements and those for the difference of 
the diagonal elements coincide with the antitrace map. 
Thus, from the trace and antitrace map, we can deter- 
mine any physical quantity related to the global transfer 
matrix of the system. Finally, as examples of applica- 
tions of the trace and antitrace map method, we inves- 
tigated the transmission problem for optical multilayers, 
harmonic chains, and electronic systems arranged accord- 
ing to the FS or the TMS. 

The trace and antitrace map method developed here 
can be expected to have many applications in the study 
of one-dimensional aperiodic systems. 
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APPENDIX A: RELATIONS FOR UNIMODULAR 
MATRICES 

For convenience, we present a collection of relevant 
identities, which are used in the construction of the trace 
and antitrace maps in Sees. [0], III, and [V. 

The nt h pow er of a unimodular 2x2 matrix A can be 
written asll3E2 

A n = U n (x A )A-U n - 1 (x A )I, (A.l) 

where / is the unit matrix and 
XI - A™ 



U n {x A ) = 



A+ = 



A-i 



A- 



X A ± yG 



(A.2) 



Here x A and A± denote the trace and the two eigenvalues 
of A, respectively, and A+A_ = Act A = 1. The functions 
U n (x) are related to the Chebyshev polynomials of the 
second kind C n (x) by U n (x) = C n -i(x/2). From the 
definition of the functions U n (x), it follows 



-1, U (x) = 0, 
1, U 2 (x)=x, 



U-i(x) 
U x {x) 
U 3 (x) 

U n +l(x) = xU n {x) - U n -l(x), 

U 2 n (x) = U n+ i(x)U n -i(x) + 1 



x 2 - 1, U 4 (x) = x 3 



2x, 



(A.3) 



In order to study the antitrace maps, we need the fol- 
lowing identityEj 



y A B = x B y A + x A y B - Vb A 



(A.4) 
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for the antitraces of two unimodular 2x2 matrices A and 
B. Now, we briefly prove this identity by introducing an 
auxiliary matrix 



1 = 





-1 



7 2 = -I, det( 7 ) = 1. (A.5) 



For the matrix A, we have 

Va = x Al . 
Then the antitrace of AB is given by 

Dab = x ABl = —x lAllB . 
Let A, B, and C be unimodular matrices. ThcnEI 

xabac = xabxac + x B c - x B x C - (A 



(A.6) 
(A.7) 



Applying the above identity to Eq. (A.7) and using 
Eq. ( |A.6| ) again, we obtain Eq. (A.4). 

It should be pointed out that Eq. (A.4) is valid for any 
pair of 2 x 2 matrices, and it follows directly from the 
identity 

AB = (x AB - x A x B )I + x A B + x B A- BA (A. 9) 

which holds for any pair of 2 x 2 matrices. The detailed 
proof of this identity can be found in Ref. Here, we 
only need to consi der unimodular matrices. 
For n = 2, Eq. (A.l) becomes 



A 2 = x A A - I. 



(A.10) 



This is the well-known Cayley-Hamilton theorem. From 
the theorem, we have 



.4 



xba - 
Vba 



+ A- 1 = x A I, 

X A -l = X A , 

xba- 1 = xaxb, 

Vba- 1 = x A y B , 



(A.ll) 



From Eqs. (A. 9) and (A.ll), we can prove the following 
useful relations 

BAB = A- x A I + x AB B = x AB B - A' 1 , 
BA~ X B = (x A x B — x AB )B — A. (A.12) 



Finally, from Eqs. (fO|), flO| ), and flATOj ), we obtain 



the following relations 



X A 2 



2. 



x A n = U n+ i{x A ) - U n -i(x A ), 
y A 2 = x A y A , 
y An = U n (x A )y A . 

This completes our collection of identities. 



(A.13) 



APPENDIX B: ANTITRACE MAPS FOR SOME 
METALLIC MEAN SEQUENCES 

The trace and antitrace maps for the golden mean and 
the copper mean sequences were discussed explicitly in 
the main part of this paper. Here, we give the trace 
and antitrace maps for some other prominent examples 
of metallic mean sequences. 

From Eqs. ( |2.26| ), Q2.27| ), and ( |Q| ), the trace and an- 
titrace maps for the silver mean case (to = 2, n = 1) are 
obtained as 

xi+i = — —[x^x 2 ^ - 1) - Xi- 2 ] - xi-i, (B.l) 

Xl-l 

xi 



yi+i 



Xl-l 



-(yi-2 + yi) + {x t - l)yi- 



(B.2) 



For the bronze mean sequence (m = 3, n = 1), we find 



xf-1 

xi+i = — - 

x i-i - 1 



[xi(xf_ 1 - 2a; i _i) - x^ 2 ] 



-XlXl-l, 



(B.3) 



yi+i = 



-{yi-2 + xi-iyt) + (xf - 2xi)yi-i. 

(BA) 



Finally, for the nickel mean case (to = 1, n = 3), the 
result reads 

Xl + l = (xf_ 1 - l)(xixi_i - xf_ 2 + 3x ; _ 2 ) 

-X/£C;_i, (B.5) 

vi+i = (xf-i - 1)(^_2 - !)y/-2 

+x l (xf_ 1 - l)yi-i - xi-iyi. (B.6) 
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TABLE I. The initial conditions for the trace and antitrace 



map. 





m = i 


ra — I 


m = j 




m — 4 


Xl 














x 2 





-0i 





0i 


V2 


-0i 





01 





yi 


2 


2 


2 


2 


V2 


01 





-0i 





W2 





-0i 





0i 



TABLE II. The trace, antitrace, and transmission coeffi- 
cients for FC(m) with m = 2q. The upper (lower) signs refer 
to even (odd) values of q. 



I xi yi tj_ 

10 2 1 

2 ±ei 4/ gl 

3 o i/em 4/ e 2 m 

4 ei o 4/el 

5 l/02m 4/gl„ 

6 ±01 4/0? 

7 l/0 3m 4/0§„ 



TABLE III. The trace, antitrace, and transmission coeffi- 
cients for FC(m) with m = 2q + 1. The upper (lower) signs 
refer to even (odd) values of q. 



I xi yi U 

10 2 1 

2 I/01 4/0? 

3 -0i 4/0? 

4 =Fl/0m+i 4/ 

0m + l 

5 l/0 m 4/0^ 

6 0i 4/0? 
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